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Abstract 

We develop a factorization method for q-Racah polynomials. It is inspired by the approach to 
q-Hahn polynomials based on the g- Johnson scheme but we do not use association scheme theory nor 
Gel'fand pairs, but only manipulation of (/-difference operators. 

El 

1 Introduction 

The factorization method is a powerful tool for investigating orthogonal polynomials. In the treatise [T5] . 
W. Miller summarized his work in which the factorization method for classical orthogonal polynomials of 
a continuous variable was developed in the framework of representation theory of Lie algebras. The same 
author treated the discrete case in [301 122] and a q-analog in [3T] . More recent papers in this direction 
are [2111] and [T7J. On the other hand, in the theory of the q- Johnson scheme, developed by P. Delsarte 
[7J[5], C. Dunkl [pj and D. Stanton [27J (see [5], chapter 8, for a recent account) a factorization method 
for q-Hahn polynomials is implicit. Indeed, in a previous work |26) we have introduced a factorization 
method for multidimensional q-Hahn polynomials in order to develop a tree method for them. In the 
one dimensional case it coincides with that in the works of Delsarte, Dunkl and Stanton but it does not 
involve any algebraic structure but only the manipulation of suitable (/-difference operators. 

The aim of the present paper is to extend the factorization method for one-dimensional q-Hahn poly- 
nomials in 26 to the q-Racah polynomials [2] [12l [14j [15l [16l [23] . This yields a new approach to these 
polynomials, based on the manipulation of a lowering and a raising q-difference operator and a natural 
factorization of the second order q-difference equation satisfied by the q-Racah polynomials. By elemen- 
tary linear algebra methods we obtain the solution of this equation and the corresponding orthogonality 
relations, with a 605 nonstandard expression for the q-Racah polynomials. Then we use a basic hyperge- 
ometric transformation formula to obtain the classical 4^3 standard expression. Another transformation 
formula is used to obtain the norm of the polynomials. At the end of the paper we show how our calcu- 
lations simplify when we take q — > 1, that is when we consider the Racah polynomials. 



The classical approaches to the q-Racah polynomials require to develop general theories (of basic 
hypergeometric series in |12j or of difference equations on non- uniform lattices in |23j). that clearly give 
more general results in other directions. Our method is more direct and elementary, because it is focused 
on a particular analysis of the structure of the finite difference equation satisfied by those polynomials. 

Following a suggestion of the anonymous referee, we begin our paper with an abstract formulation 
of our version of the factorization method (Section [5]) and then we show the route to follow to develop 
factorization methods for Hahn, q-Hahn and Racah polynomials (Section [3]). We point out again that 
the theory in Section 2 apply also to the multidimensional case [26] (but in this case also a tree method 
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must be developed). 

It is intriguing to know if the multidimensional Racah operators introduced by J.S. Geronimo and P. 
Iliev in [13] admit a factorization method that generalizes our one-dimensional case (for q=l), along the 
line of the factorization method for the multidimensional g-Hahn operator in [26 (for q = 1, see also [2"5]V 
In the case of positive answer one might investigate for a g-analog. It should also be interesting to derive 
the factorization method in the present paper from the fact that the g-Racah polynomials are connection 
coefficients for different families of bidimcnsional q-Hahn polynomials: see [11] and our generalizations in 
[25] . The calculations in Section 6 of [53] (in the q = 1 setting) suggest that this should be possible. 



2 The factorization method in an abstract form 

In this Section we present our version of the factorization method in an abstract form. A similar but 
slightly different algebraic formulation is in [5D]. Suppose that Vo, Vi, V2, ■ ■ ■ , Vjv> • • • is a sequence of 
finite dimensional complex vector spaces, with dim Vn = N + 1. Suppose also that Vjv is endowed with 
an hermitian scalar product (•, -)v N - Let 

T>n : Vn — > Vn 
be a linear selfadjoint operator, N — 0,1,2,.... 

Definition 2.1. A factorization method for the sequence of operators T>n, N = 0,1,2, ... is given by 
two sequences of linear operators 

K N ■ V N -> V N +i, N = 0,1,2,... and C N : V N -> V N -i, N = 1,2,3... 
(raising and lowering operators) such that: 

{K N - 1 h,h)v N = {h,£Nh)v N - 1 , forallf 1 €V N - 1 ,f 2 eV N ,N>l, (1) 
that is Cn is the adjoint ofR-N-\ } 

n N -iCN = -VN + a N lN, N>1, (2) 

and 

Cn+i'R-n = —T^n + onIn, N>0, (3) 
where aN,bN are constants and In the identity operator on Vn- 



From ^ and ([3]), we get the fundamental commutation relation: 



Cn+iR-n — TZn-i£n — {on ^ Q,n)I_ 



N ■ 



N > 1. (4) 

Note also that the selfadjointness of T>n is also a consequence of |T|) and ©. It is convenient to set 
H-i = 0, Cq = and define do in such a way that T> = a I, so that ^ and (0| hold true also for N = 0. 
Now we give some basic consequences of these identities. 

Lemma 2.2. For N > n > 0, we have: 



CnTZ-n—iR* 



N-2 ■ ■ 



■n n = n 



Af-2^JV-3 



'JV-1 



^2(h - Ofc) 



_k—n 



IZN-2HN-3 ■ 



■n n . 



(5) 
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Proof. For n = N — 1 it coincides with Q (with TV replaced by N — 1). The general case can be proved 
by induction on TV — n: 

Cn'R-N-i'R-N-2 ■ ■ ■ Tin = ^JV-2-^JV-l^iV-2*^JV-3 ' ' ' Tin + ipN-X — a TV -l) 7?- TV -2 7^ TV -3 ' ' ' Tin 

[N-l 



TIn-2TIn~3 ■ ■ ■ Tln-lCr, 



TZn-2TIn-3 ■ ■ ■ Tin, 



where the first equality follows from (Q| and the second from the inductive hypothesis. 
Lemma 2.3. IfO<n<m<N and C n f — 0, then 

'N-l h 



Cm+lCm+2 • ■ ■ C-nHn-iTIn-2 1 ' ' Tlnf — 



_h—m k—n 



m— 1 >^m — 2 ' ' ' 



Tlnf. 



□ 



(6) 



Proof. For m = N — 1 it is an immediate consequence of ([5]) ; the general case follows by induction on 
N -to: if £„/ = then 



£-m+l£m+2 • • • £nTIn-iTIn-2 ' ' ' Tlnf — 



N-l 



_ k—n 

N-l h 



m + 1 £ m + 2 Cn-iTIn-2TIn-2 ■ ■ ■ Tl n f 



n E(^-^) 



Tlm-lTlm-2 ' ' ' 



_h—m k—n 

where the first identity follows from ([S]) and the second identity from the inductive hypothesis. 

Lemma 2.4. IfO<n<m<N, f\ € Ker£„ and /2 € Ker£ m , i/ien 

nv-i ft 



□ 



n e^-^) 



(/i,/a>i 



{TIn-iTIn-2 ■ ■ - Tlrnfl,TlN -{TIN -2 ' ' ■Tlnf2)v N — 5 n 

_h—n k—n 

Proof. It is an immediate consequence of (QJ and (|6]), noting also that (p} implies that Ran7?. m _i_LKer£ 

Corollary 2.5. Suppose that N > n. Then for all f € Ker£ n we have 

pv-i ft 



(7) 



□ 



\Hn-iTIn-2 ■ ■ ■ Tlnf\\v N — 



n E( 6 *-o 



,h=n k—n 



\ffv n 



(8) 



In particular, the linear map 



Ker£„ — ► F w 

/ 1 — > TIn-iTIn~2 ■ ■ ■ Tlnf 



is infective. 

Lemma 2.6. If N > n > and f £ Ker£ n , then 



T>N [Hn-iTIn-2 ■ ■ ■ Tlnf] — 



N-l 

a N - 2J ~ Ofc) 



TIn-iTIn-2 ■ ■ ■ Tlnf- 



(9) 



Proof. Indeed, we have 



AT-l 

aN - - a k ) 

k—n 



n N -{R 



N-2 ■ ■ ■ 



where the first equality follows from and the second equality from © (or |(SJ)). 



□ 



Theorem 2.7. Suppose that Ker£„ is non trivial for n > and choose ip n € Ker£ n; tp n 0. Set 
fN,n = TL N -i1l N -2 ■ ■ ■ UnVn ■ Then 

fN,0, /jV,l) • ■ • i /iV,n) • • • ) /jV,JV-l) /jV.JV 
is a complete set of mutually orthogonal eigenvectors for T>n , N > 0. Moreover we have: 



rw-i 



/jv-l,nj 



(10a) 
(10b) 



Proof. Indeed, (O ensures that each / n is an eigenvector of ©at, while ([7]) ensures that these eigenvectors 
are mutually orthogonal. Recall also that dim Vn — N+l by hypothesis. Moreover, (|10a[) is a consequence 
of the definition of f^, n while (jlObj) is a consequence of (J6|) (or ([5|)). □ 

3 Factorization of second order finite difference operators 

If we want to apply the factorization method developed in the previous Section to the classical families of 
orthogonal polynomials, the main problem is to find sufficient conditions for the factorization of a second 
order finite difference operator. In this Section we give sufficient conditions for the existence of first order 
difference operators that satisfy the relations and ([3]) ; these conditions apply to Hahn, q-Hahn, Racah 
and g-Racah polynomials; the condition ([1]) will be verified directly in the case of g-Racah polynomials 
in Proposition 15.21 

Suppose that V N is the space of all functions / : {0, 1, 2, . . . , N} — > C (with /(-l) = f(N + 1) = 0) 
and that T>n is a second order difference operator of the form: 

V N f(x) = B N {x)f{x + 1) - [B N (x) + D N (x))f(x) + D N (x)f(x - 1), 
where Bn,Dn are defined for x = 0, 1, 2, . . . , N. 

Proposition 3.1. Suppose that Bjy and Dpj may be factorized in the form 

B N (x) = B ( V (x) ■ B% } (x) D N (x) = Z?£> (x) ■ D%> (x) 

where B$ , Bjp , ( defined also for x = — 1) and D^- ( defined also for x = N + 1) satisfy the relations: 



B%l 1 (x)B$_ 1 (x-l)=B N (x), 
D^ ) _ 1 (x~l)D^_ 1 (x)=D N (x), 

Bfil^x - l)D$_i(aO + D^l^B^ix - 1) - B N (x) - D N (x) = a N , 



D^'{x)B^ ) {x-l)+B^>{x)D^'{x + l) - B N {x) - D N {x) = b N , 



l(2 



(11a) 
(lib) 
(11c) 
(lid) 
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for x = 0, 1, 2, ... , N, and a^, on are constants. If we define first order difference operators by setting 



K N f(x) = -B%\x - l)f(x) + Df(x)f{x - 1), 
C N f(x) = -D^l^fix) + B^Wfix + 1), 

then IZn arid Cn satisfy the conditions ^ and (J3J (with ajv and 6jy given by (| 1 lc[) and f| 1 ld[) ). 
Proof. We have: 



n N .xL N f{x) = - Dg^(i - VDfil^xMx - 1) 



B^x VDfil^x) + DXL^B^x - l)j f(x) 
= - T>xf(x) + a N f(x), 

where in the first identity we have just applied the definition of 1Z n and Cn, and then we have used 
(|TTal) . (|llbl) and ([TTc]) . The proof of © is similar. 

□ 

It is convenient to set T)n = £n = 0. 



(i) 



>(2) 



If the functions fN,n hi Theorem 12.71 are constructed in the setting of Proposition 13. 1[ then the 
formulas (1101) take the form of first order difference relations: 



(2) (2) 
B N \ X ~ 1 )InA x ) + D N ( x )fN,n(x - 1) = fw+l.n, 



jV-1 



a-kj 



fN-l,n(x), 



(12a) 
(12b) 



Now we indicate factorizations techniques for Hahn, g-Hahn and Racah polynomials. The details for 
Hahn and g-Hahn are straightforward and are purely analytical versions of the classical approaches of 
Delsarte [8] , Dunkl [10] and Stanton [27] ; an exposition of the results in these papers may be found in 
[BJ, Sections 6 and 8. The g-Hahn case may be found in [2B] (where the method is applied also in the 
multidimensional case). We sketch some details for the Racah polynomials and in Section [7] we show how 
the q- Racah case simplify when we make q — > 1. 

Hahn polynomials. 

In the case of Hahn polynomials Q] [T4] [X5l [SI [23] we have: 

B N (x) = (x + a + l){x-N) and D N (x) = x(x - f3 - N - 1). 

If we take 



B N \x) 



■a+1, 



B%\x)=x-N, 



P-N -1, 



D${x) 



it is straightforward to check that the relations (fTTj) are verified, with etjv = N(a + (3 + N + 1) and 
b N = (N+l)(a+*3 + N + 2). 
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q-Hahn polynomials. 

In the case of g-Hahn polynomials [21 [TSl EH [23J we have: 

B N (x) = (l-aq x+1 )(l~q x - N ) and D N (x) = q^-^l - q x )a(Pq N+2 - q x+1 ). 
If we take 

B^(x) = l-aq x +\ B^(x) = l-q x - N , 
D$(x) = a^* 2 - <f = ^(l - q% 

it is again straightforward to check that the relations (fTTj) are verified, with ajy = — (1 — a^g +1 )(1 — q ) 
axidb N = -(l-q- N - 1 )(l-a/3q N+2 ). 

Racah polynomials. 

In the case of Racah polynomials [T^l [TSJ [23] we have: 

(x + a + l)(x + /3 + S+ l)(x - N)(x + S - N) 



B N {x) = 
D N (x) 



(2x + 8 - N)(2x + 8 - N + 1) 
O-a + 5- AT - l)(x — f) — N — l)x{x + 8) 
(2x + 8 - N)(2x + 8 - N - 1) ' 



We can take 



(i) _ (a: + a + l)(a: + £ + <? + !) (2) _ (x-N)(x + S-N) 

** [x) - 2x + 8~N ' N [X) ~ 2x + 8 - N + 1 : 



(1) = (x-a + J-7V-l)(x-/3-jV-l) (2) = x(x + £) 
^ V ; 2x + 5 - TV ' N K } 2x + 5-N-l' 

We just sketch the proof of (|1 lci() : 



D%Xx)B%\x - 1) + B%\x)dP(x + 1) - Bjv(^) - £>jv(z) 
(x-a + S-N- l)(x - - N - 1) 



(2x + <5-iV)(2a; + 5-iV- 1) 
(.T + Q + l)(a: + /3 + (5 + l) 
+ (2x + (5-iV)(2a; + (5- AT + 1) 
7V + l 



[(x -N- l)(x + 8-N-l)-x(x + 8)] 
[(x + l)(x + 8 + l)-(x- N)(x + 8-N)] 



(x - a + 8 - N - 1)0 - P - N - 1) + + a+ l)(x + + 8 + 1)] 



2x + 8 - N 1 

= {N + l)(a + P + N + 2). 

Therefore 6jy = + l)(a + /3 + 7V + 2). Similarly one can check the other identities (fTTj) and that 
ajv = 7V(a + ^ + iV + 1). 

4 The difference equation of g-Racah polynomials 

Let a, /3, 7, 5 be complex parameters, fix a real number < q < 1 and set 

(1 - aq x+1 )(l - P8q x+1 )(l - iq x+1 ){l - -ySq x+1 ) 

(1 -7<5 g 2 ^+ 1 )(l -j8q 2x + 2 ) 
g(l - g x )(a - 7 8q x )(P - lq x ){l - 8q x ) 




G 



Let JVbea positive integer. The q-Racah operator T> = T> a ,p,-y,s is the second order g-difference operator 
defined by setting 

Vf{x) = B(x)f(x + 1) - [B{x) + D(x)]f(x) + D(x)f(x - 1) 

for each complex valued function / defined for x 6 {0, 1, . . . , iV}. The g-Racah polynomials are the 
eigenfunctions of T> when one of the following conditions is satisfied: 

— N as — jV -N 

aq = q , poq = q or 75 = q 



Consider the following substitutions on the parameters a, /3, 7, S: 

o-(a,/3, 7 ,8) = (a,p,p6l) and f(a,ft 7 ^) = ( 7l ^,a^). (13) 

p 7 a 

It is easy to see that T> is invariant under a and r. Moreover, a switches the condition f35q = q~ N with 
7<7 = q~ N and fixes aq = q~ , while r switches aq = q~ with 75 = q~ N and fixes f3Sq — q~ N . Note 
also that, if we think of a and r as transformations acting on C 4 \ {(a, /3, 7, 5) : af3j6 7^ 0}, then tr and 
r generate a group isomorphic to the symmetric group on three objects (a 2 = r 2 = (err) 3 = 1). In what 
follows, we will concentrate on the case 7*7 = q~ N . Indeed, using the symmetries (I13p . our results can be 
easily translated into the cases aq = q~ N and (35q — q~ N . 

The space Vm of all complex valued functions defined on {0,1,..., N} will be endowed with the 
following scalar product: 



if f \ _ N(N+i)/2 { a6 1 q,Pq;q)N ^ (aq,(35q,q N ,Sq N ;q) x (1 - 5q 2x N ) —— 

{h ' h)VN q (q,6-i;q) N ^ (q,a^6^ , /3"V JV ,Sq;q) w ' (a(3q) x (l - 5q-^) h[X)h[X) - 

(14) 

It is the usual scalar product for the g-Racah polynomials in the case jq = q N (see [TH HH IT5"]). 
multiplied by g jV ( jv + 1 )/ 2 ("^ , An equivalent but useful expression is the following: 

(fuh)v N = £ (^"^"^"'ig)^^- 1 ^)^-- . r ^(^i)/ 2+ , (1 _ 5-\ N -**)h{x)Mx). (15) 
£r£ {q;q)x{q;q)N-x{8 1 q x ;q)N+i 



The expression (|15l) may be easily obtained from (| 14[) by m6a/ii of the elementary identities -r ■■■ ]■-■■ -x~ , t , 
(_l) V ^(JV-x)+,( :c+ i)/2 (ag . g)iv _ :c) (5-i ;9f ) w (i_^-W) = -^(^^^(l-J-i) and (lg;g)I 

a x g :c ( :t: + 1 )/ 2 (-l) x (a-V" :c ;5)x. 



5 Factorization of the g-Racah operator 

For jq — q~ N we have: 

(1 - aq x+1 )(l ~ f38(f +1 )(l - q x - N )(l - 5q x - N ) 



B N (x) 
D N {x) = 



(l-5q 2x - N )(l-6q 2x - N + 1 ) 
{a - 6q x - N - 1 )(f3q N + 2 - q x+1 )(l - q x )(l - Sq x )q- N - 1 



{l-Sq 2x - N )(l-Sq 2x - N - 1 ) 
In the notation of Proposition 13. 1[ we may set: 
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(1) _ (2), _ (l- q x-»)(l-Sq*- N ) 

°N W ~ 1 _ § q 2x-N ' W — 1 _ fa2x-N+l ' 

(1) _ (a-Sq*- N -l)(l3 q N +2- q ^) (2) _ (l-g^g-g^g-^- 1 
The corresponding raising and lowering operators have the form: 



ftjv/fr) = 1 _ ^zjvri /(a) + 1 _ Sq 2x-N-i /( x - !)« 

( a _ Sa x ~ N )(Ba N+1 - n x+1 ) (1 - aa x+1 )(l - 85a x+1 ) 



(16) 



for all / G V/v- Observe that if / G V/v is constant, f = c, then £/v/ = (1 — a/3g A ' +1 )c, 72. jv/ 
-(1 - q-t-^c and XV/ ee 0. 

Proposition 5.1. TTie functions defined above satisfy the conditions (|lip. witt ojv = —(1 — q~ N )(l 
a(3 q N+1 ) and b N = -(1 - - aPq N+2 ). 

Proof. For instance, we prove (II lc|) . From the elementary identities 

(a - Stf-^tfq 1 ** 1 - - (1 - aq x+1 )(l - f35 q x+1 ) = -(1 - Sq 2x - N+1 )(l - a/3q N+1 ) 

and 

(1 - a q x )(l - /3S q x ) -{a- < 5< 7 a; - Ar - 1 )(/3g Ar + 1 - q x ) = (1 - Sq 2 ^' 1 )^ - af3q N+1 ). 

we get: 

B^L^x - ^D^l^x) + D^l^B^ix - 1) - B N (x) - D N (x) 

(a ~ 5 q x - N )(j3 q N+1 - q x+1 )(l - q x - N )(l - S q x - N ) 
- ' B N (x) 



{l-5c? x - N + l ){l-5 q 2x - N ) 
(1 - aq x ){\ - j35a x ) q - N (l - q x ){\ - 5q x ) 
+ (1 - Sq 2 *-"- 1 )^ - 8q 2 *- N ) 

1 - a(3 q N+1 



D N {x) 



I ( - ,, y [(1 - q X )q- N {^ ~ Sq x ) - (1 - q x - N )(l - 6q x - N ))] = -(1 - q- N )(l - aBq N + l ). 

In the same way one can prove (|lld[) and get the value of 6jv- d 
Now we check ((T|). 

Proposition 5.2. For f\ e V/v-1,/2 € V/v, -TV > 1, we fterae: 

(R.N-ifi, h)v N = (fi,Cffh)vs-n (17) 

i/iai is £at zs the adjoint ofR-N-i- 

Proof. Using the expression (|15[) for the scalar product, after some immediate simplifications, we have: 
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JV-1 



I'D t f\ (gglft _ g ^)x(/3g,aS 1 g;g) N -x rx „n(n+3)/2 

{i<-n-iJi, h)v N - > 7 7 s tfzi -x — \ '" q 

^ {q;q)x{q;q)N-x(d l q x ;q)N+i 



x=0 



x S-\l - q x - N ){\ - 5q x - N )h(x)h{x) 
(m, P~ 1 8~ 1 q~ x ;q)x(l3q, aS- 1 q; q) N ^ x _ ^ XqN ( N+1) / 2 ^ 3 



x&-\l-<f){l-6<f)h(x-l)h(x). 

Similarly 



It r f \ \ " ^ g ;q)x(Pq,aS 1 q;q) N -x-l RX m(N+1)/2-x-1 v 

^ wgjxtegjjv-x-x^ q ,q)N 

X ^ X (Q - 6q x - N )(f3q N+1 - q x+1 )h{x)Ux) 
_ (<*q, P~ l 5~ l q~ X iq)x{Pq, ad^q; q) N -x-i _ ™jf( W )/ 2 -»-i x 

x ^(l - a<f +1 )(1 - ft I+1 )/i(i)/ 2 (i + l). 
Then we have (R. N -ih, f2)v N = (fi, C N f 2 )v N - 1 because 

(P^aS-^N-s-i JV(JV+1)/2 _ Z _ 1 _ 5(f -iv )(/? iv+i _ x+i) 

= (Pq,a6- 1 q;q) N - x _ N _ x)( c-i - X+N) N(N+i)/ 2 -x-i (a - Sq x ~ N )(/3q N+1 - q x+1 ) 

(q;q)NM5- 1 q- x ;q)N+i [ q )[ q )q (1 - a6^q N -)(l - j3q N -) 



((3q : aS 1 q]q) N - 

(q;q)N-x(5- 1 q- x ;q)N+i 



-x q N(N+3)/2- x ^ _ qX-N^ _ Sq x- N) 

wvjn-xk." -q ~,q>N+i 

and 



(aq,f3 M 1 q x ;q) x (f3q,aS 1 q; g) N -x-i ^ Xq N(N+i)/2- x -i^ _ at f+lyi _ (3Sq x+1 ) 

(<?; q)x(q; q)N-x-i(5- 1 q- x ;q)N 



x=y-l 



(q;q)y-i{q;q)N-y{5 L q y+L ;q)N 

= (gg) /g'M- 1 ^-^; gp^Cgg, a^" 1 ^; g)jy- y (1 - - ^iV+lVa-j//-, _ jn 

(q;q)y(q;q)N- y (S- 1 q-y;q) N+ i ' 1 - V ' P ? ^ PQ> 

= {aq p-^q-y-q) Wq aS^q;q) N ^ y py qN(N+1)/2 - y 

{q;q)y{q;q)N-y{o L q v ;q)N+i 

□ 

From the Propositions above we can deduce that the sequences of operators TZn and Cn in (I16p form 
a factorization method for the g-Racah operators (Definition I2.1j) . 

Lemma 5.3. The operators T>m, "R>n and Ln of this Section satisfy all the properties proved in Section 
2 in the abstract setting, with 
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N-l 



N-l h 



k—n 

II E( 6 * - a *) = 9- (JV - m)(JV+m+1)/2 (<z m -" +1 ;<z)A r - m Mg n+m+2 ;g) 

-n 

V-l /i 

II - a *) = f (N -" p+ " +1)/2 (9;9),-»W? 2 " +2 ;9)«-„ 



N—mi 



h—ra k 



h—n k—n 



N-l 



/c— n 

Proof. These are elementary calculations. To simplify just note that a& + i = 



(18) 

(19) 
□ 



We end this Section giving an explicit form for the repeated action of the operators TZn- In what 
follows, as in [9l [TTj . for a, 6 € K we set a V 6 = max{a, 6} and a A 6 = min{a, 6}. We also recall the 
definition of the Gaussian or g-binomial coefficients: for integers k, n, with < k < n, 



(q n - k+1 ;q)k 



q («;?)* 

Lemma 5.4. Suppose that < n < N . Then for all f E V n and < x < N we have 



1Z N -iTl N -2---'R- n f(x) = (q;q) N - n E 

y=(x-N+n)V0 



' N 


— X 




X 


n 


- y _ 


Q 


y. 



(Sq x - N ;q) v (5q x - n +y +1 ;q) n . 



(6qv- n ;q) n+1 

x (1 - 5q 2v - n )q^ y+N - x -^q-^ N - n ^ N+n+1 ^ 2 f(y). (20) 

Proof. The proof is by induction on N. For N = n it is trivial; it is also easy to check that for N = n+1 
it coincides with the definition of 7Z n . Suppose that (|20[) is true. Then, with a repeated use of the 

m 
k 



m + 1 
fc 



elementary identities 

ll N K N -i ■ ■ ■ TZ n f(x) = (q; q) N . 



i-g" 



and (a;q)fc+i = (a; q)k(l - aq k ), one can get: 



^ E 

t,= (a;-Af+n)V0 



X (1 - ^2 2 /-n\ ?! /(j/+iV- a :-n) ? -(iV-n)(Ar+n+l)/2 . 

(x— l)An 
+ (q;q) N - n E 

y=(x— N+n— 1)V0 



"JV- 


a; 




a; 




n — 


y . 




y. 




-l)/2 . 


5 a 


-jv-i (1 



(Sq - N ;q) y (Sq x - n +y+ 1 ] q) n - l 
(Sqy- n ;q) n +i 

- x+N+1 )(l-Sq x - N - 1 ) 



~N -x + l~ 




X-l 


n-y 


Q 


y 



1 - 5q 2x ~ N - 1 
(dq x - N - 1 :q) y (Sq x - n +y-q) n _ l 



f(y) 



{8qy- n -q) n+1 



(1 - Sq 2 y- n )q v(v+N - x - n+1) q- (N - n)(N+n+1)/2 q- N - 1 ^ f(y) 



(q;q)N- n +i E 

y=(x— N+n— 1)V0 
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i - <v 

(^- Ar - 1 ; g ),(^-W 1 ;< Z ) n _, 



(*g»-";<f) 



n+l 



X (1 - § q 2y-n\ q y(y+N-x-n+l) q -(N-n+l)(N+n+2)/2 x 
(1 - gJV-x-n+y+l^j _ fiqX+y-N-ljqX-y + (1 _ ^-^(l - 5^+2/-") 



/(</), 
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and the last expression coincides with ([20)1 for IZ^TZn-x • • ■ 1Z n f(x), because the fraction in square 
brackets is equal to 1. 

□ 

6 The g-Racah polynomials 

Now we are ready to derive the q-Racah polynomials r n (x; a, (3, S, TV; q) as an explicit complete orthogonal 
system of eigenfunctions for the g-Racah difference operator (Theorem |2.7|l . We will first give a non 
standard 60s expression for r n (x; a, (3, 5, N; q). Our machinery reduces the computation of the norms of 
these polynomials to the case N = n, which requires the following summation formula for a terminating 
very-well-poised series (formula (2.4.2) in |12|): 



a,qy/a,—qy/a,b,c,q " aq 
/a,—^/a,aq/b,aq/c,aq n+1 be 



n+l 



(aq,aq/bc;q) n 
{aq/b,aq/c;q) r 



(21) 



In order to get the standard expression for the q-Racah polynomials, it is convenient to use the normalized 

map ^ q ("-nw+n+Wn N _ 1 n N _ 2 ---n n . 

Theorem 6.1. For < n < N set 

r n (x; a, 0, S, n; g) = {a[3q n+l T ^ /'""o'T^P q ^ , x = 0,l,...,n, (22) 

(aq,f3dq;q) x 



r n (x; a, (3, 6, N; q) = Ml^-XWIJ/^k^^^^ . . . H n r n (x; a, (3, 6, n; q), x = 0, 1, . . . , N, 



Then we have the following explicit formula 



r n (x;a,(3,S,N;q) = 



(q N - n+1 ;q)r 



E 



y=(x-N+n)\/0 



' N — 


X 




X 


n — 


y _ 


Q 


y . 



(6q*- N ;q) y (5q 



x— n+y+l . 



q)n- 



(23) 



{8qy- n -q) n+l 
„—lX„-n a—l„-n. 



(1 - 5q 2 y- n ) q v(y+ N -*- n ) x (apq n+1 ) y ^ "'^ 4 — ■ (24) 

{aq,f35q;q)y 



Moreover, the functions r n (-;ct,/3,S,N;q) satisfies the following q-difference equation: 

V N r n {x- a, (3, 5, N; q) = gT n (l - q n ){l - a[3q n+1 )r n (x; a, /?, S, N; q), 
and the following orthogonality relations: 



(25) 



(r n (-,a, (3, 5, N; q), r m (-, a, f3, 6, N; q))v N 

- A „N(N+i)/2 rx „-N^n (a(3q n+1 \ gWg, aS^q, (3q; q) n 

{ q ' {l-a[3q^){q ]q ) N { q -»,[36q,*q-, q ) n {2b > 

Proof. First of all, note that r n (x; a, (3, S, n; q) (as defined in (|22p) is the unique function in V n satisfying 
the first order difference equation C n (p — 0, with the initial condition ip(0) — 1; indeed, £ n tp = is 
equivalent to 
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(a - Sq x - n )(pq n+L - q x+L )<p(x) = (1 - aq x+L )(l - f35q x+L )ip(x + 1), x = 0, 1, . . . , n - 1. 

Then (pH|) immediatly follows from (|20[) . while (|25l) follows from © and (|19p. The orthogonality relations 
(f2"B]) follows from the selfadjointness of 2?jv; to compute the norm, we first derive the case N — n. Using 
the expression (fT^|) for the scalar product, the explicit formula (|22p and the identity 

l-Sq 2x - n _ (5q- n+2 ;q 2 ) x _ (VSq- n / 2+1 , -y/6q~ n / 2+1 ; q) x 



1-Sq~" {5q- n ;q 2 ) x {VSq- n / 2 ,-VSq- n / 2 ;q) x 
one can easily get the following 6</>5 expression for the norm of r„(-; a, (3, 6, n; q): 

||r„(-; a, 0, S, n; q)\\ 2 Vn =g » (w+1}/a \ q)n x 



X 605 



W"/ 2 , - W l/2 , /3<5g, ag, <5g ' 9 ' ap<Z 



Applying (f2"Tj) with a = <5g n , b — /3 1 q n and c = a 1 <5g ™, we can sum this e4>5 expression, obtaining 

r n (-;a,^,d,n;g) f =g 1 + " -— r 

[Q,S 1 ;q)„ ((35q,aq;q) n 

=q n/_ S \n ^ a/?g n+1 ; g) n 



Finally, using ((SJ) , (TT51) and ([2"31 (and the elementary transformation formulas < - 9 ' g - )jv " — — — 



(q,f35q,aq;q) n 

M-wm 1 o o . . , 



and (a/3g™ +1 ; q) n (a/3q 2n+2 ; g)jv-« = ^"f_ a/ 3 g ^»-fi +1 )i one can S e t easily the expression for the norm of 
r n (-;a,P,6,N;q) in (29 . 

□ 



It is easy to check that the expression for the norm of the g-Racah polynomials in (|26p coincides with 
the usual expression (see [13 El [15]) multiplied by the factor q N(N+1),2 ^j^r^ L - 

Corollary 6.2. r n (x; a, ft, 8, N; g), n = 0, 1, 2, . . . , N, is a complete set of eigenf unctions for the operator 
V N . 



Now we indicate the action of TZn and £jv on the g-Racah polynomials. 
Corollary 6.3. The q-Racah polynomials satisfy the following first order difference relations: 

(1 - g*)(l - 6q x )r n (x -I; a, (3,6, N; q) - (1 - g^" 1 )^ 1 - Sq x )r n (x; a, /3, 6, N; q) 

= (1 - ^-^Xl - q N+1 )r n {x- a, 0, S, N + 1; q), (27) 

(1 - aq x+1 )(l - (35q x+1 )(l - q N )r n (x + 1; a, /3, 5, N; q) 

- (a - Sq x - N )(Pq N+1 - q x+1 )(l - q N )r n (x; a, /?, S, N; q) (28) 
= (1 - q N - n ){\ - apq N+n+1 )(l - 5q 2x - N+1 )r n {x- a, /?, S, N - 1; g), 
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Proof. These are immediate consequences of the relations (|12p . taking 



fa* 0=) = ( g;g )„ g (^-»)(^+"+i)/2 r " «> A «) 



(see P5J)). 



□ 



Remark 6.4. The relation ([27)1 coincides with (19) in [5] (with 7 = q~ N ~ 2 ), while (|28p coincides with 
(23) in [5] (with 7 = q~ N ). We thank A. Borodin for suggesting us that the relations in Lemma 7.3 of 
[5] could be deduced in the framework of our approach. 



In order to give the classical 4^3 hypergeometric expression for r n , we will use the following identity 
of Watson (see [12], formula (2.5.1)): 



a,q^/a,-q^/a,b,c,d,e,q ™ 
--Ja^aqjb^aqj 'c,aq/ 'd,aq/ 'e,aq n+1 ' ' 



aV+ 2 



6cde 
(aq,aq/de;q) n 
(aq/d,aq/e;q) n 



4«:S 



(7 " ,d,e,aq/bc 
aq/b, aq/c, deq~ n / a 



;q,q 



(29) 



Theorem 6.5. The q-Racah polynomials have the following 4^3 hypergeometric expression: 

r n (x;a,f3,5,N;q) = 4 <fe 
Proof. Consider the following elementary transformation formulas: 

{q N- X -n + l. q)n{q -n. q)y 



q- n ,Sq x - N ,q- x ,a(3q n+1 

aq,/36q,q- N 5 q ' q 



(9; q)r 



N-x 
n-y 



f_l\y qy{y+i)/ 2 +y(n-y) . 



{q N- X -n+l. q)i 



^l^y q y(x-y)+v(y+i)/2 . (g ^g)?/ 



x— n+y+1. 



9)n-s 



(Sq 



X — 71+1 



;<?)r 



1 - 5g 2 ^-" 



(gg-_j g),( W"/ 2+1 ; qU-VSq-^ 1 ; q) y 



(Sqv- n ;q) n+1 (Sq-^;q) n (Sq; q) y (VSq- n / 2 ; q) y (-VSq- n / 2 ; q) y 

Applying those formulas to ([24]) . after several elementary calculations one can get the following 
expression for the g-Racah polynomials: 



( N-x-n+l c x-n+l. \ 

r n (x;a,(3,S,N;q) = [Q „ . ^ ,' gj " x 



N-n+l X n -n+l 



5q- n+1 ;q))r 



X 8 7 



Sq~ n , VSq- n ^ +1 ,-VSq- n/2+1 ,f3- 1 q- n ,a- 1 5q- n ,6q x - N ,q- x ,q' n 
VSq- n / 2 ,-VSq- n ^,l3dq, aq, q N-x-n+l ^ q -n+x+i ^ Sq 



q,af3q 



N+n+2 



Setting a = Sq ™, b — f3 1 q ™, c = a 1 5q n , d = Sq x N and e — q x in the Watson identity (l2l)l) . one 
can get immediately the 4(^3 expression. 

□ 

Remark 6.6. We recall that since (q~ x ,5q x ~ N ;q) = Ylh=o(^ ~ M^)?' 1 + 5q~ N+2h ), where fi(x) — 
q~ x + 5q x ~ N , the function r n is actually a polynomial in fi(x) (of degree n); its standard notation is 

iin(//(z);a,/3,7,%) mmng. 
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7 The Racah polynomials 



Now we sketch how the calculations in Sections [5] and [6] simplify when we take q 1, that is how our 
method works for the Racah polynomials (see also the end of Section [3]) . The space Vjv will be endowed 
with the following scalar product: 



lf r \ _ (g + l)jv(a - 8 + l) N f (a + + g + 1) X (-N) X (6 - N) x ((5 -N + 2)/2) x — - 
{n ' f2)v » Nl(-S) N ^ (- a + 6 -N) x (-p-NU(5-N)/2U6 + l) x x\ 



It is the usual scalar product for the Racah polynomials in the case 7 = — N — 1 (see [Hill]), multiplied 

l) N (a-S+ 
N\(-S) N 



by ^ +1 ^T/" (5 / +1 ' >iv ■ The raising and lowering operators have the form (see Section [3} : 



(x — N — l)(x + S-N-l) ff x(x + S) 

UnI{x) = 2X + S-N-1 f{x) + 2x + 6-N~l f{x 1} 



and 



£Ar/(2;) = 2.x + ^-jV + l + 2x + S-N + l f{x + 1} ' 

for all / G Vjv Long but easy calculations show that Cn is the adjoint of 1Zn-i- In the identities 
in Section [2] now we have: Y,k=n( b k - Ofc) = (N - n)(a + (3 + N + n + 1), II^Tm Sfe=n( fe fe - a fe) = 
(m - n + l)jv-m(a + /3 + n + m + 2)Ar_ m , 

JV-1 ft 

II = (N-n)\{a + (3 + 2n + 2) N _ n , (30) 

JV-1 

a-N - - on,) = n(a + /3 + n + 1). (31) 

Lemma 7.1. Suppose that < n < N . Then for all f 6 and < 1 < # tee /iaue 



n N - 1 n N - 2 ---Tl n f(x) = (N-n)\x 

'N - x\ fx\ (x + S - N)y- x+ N- n (y + 5 + l)x- v 



X An 

X 

y=(x— JV+n)V0 



— " V n-y J \yj (x + y + 6 - N) y - x+N - n {2y + 6 -n + l) x - y 



Proof. The proof is by induction on N . For N = n + 1 it is just the definition of TZ n . Suppose that (f3"2"]l 
is true. Then we have: 
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n N n N -i---n n f(x) = (N-n)\ Y, 



N-x\fx\ (x + 5-N) y ^ x+N - n (y + 6 + l) 



x-y 



(x-l)An 

+ (N-n)\ 



y=(x-N+n-l)\/Q 



(N-n+1)) Y 



y=(x~N+n)\/0 

(N -x + l)(x + S-N -1) 
X (2x + S-N-l) 



N-x + l\fx-l 



n-y ) \yj (x + y + 5 - N) y - x+N - n (2y + 5 -n + l) x _ s 



f(y) 



(x + S - N - l) y _ x+N _ n _i(y + 5 + l)x- y -i 
n-y J\ y ) (x + y + 5-N-l)y- x+N - n -i{2y + $-n + l) x -y-- l 

xjx + S) 
X (2x + S-N-l) ny) 

N - X + l\ f x\ (x + 5 - N - l) v -x+N-n+l{y + 8 + l) x -y 



n-y ) \yj {x + y + 5- N- l) y _ x+N ^ n+1 (2y + 5 -n + l) x - y 



y=(x—N+n—l)VQ 

~(N-x-n + y+ l)(x + y + S-N-l) + (x- y)(x + y + S-n) 



(N-n + l)(2x + S-N-l) 

and this coincides with ([52)1 for IZnTZ-n-i • ■ ■TZnf(%), because the expression in square brackets is just 
equal to 1. 

□ 

In the following Theorem we will use the identity (Corollary 3.4.3 in [I]) 

a,a/2 + l,c,d,-n \ ( a + l) n (a - c - d + l) n 



5F4 • 1 

y a/2, a — c + 1, a — d + 1, a + n + 1 ' 



Theorem 7.2. For < n < N set 



r n (x;a,/3,S,n) 



(-a + 6 - n) x {-j3 - n) x 
(a + l) x (/3 + 6 + l) x ' 



(a - c + l) n (a - rf + l) r 



x = 0, 1, . . . , n, 



and 



r n (x; a, /3, 5, iV) = • • ■ TZ n r n (x; a,j3,S,n), x = 0, 1, 2, ... , AT. 



(33) 



(34) 



(35) 



TTien we /lave i/ie explicit formula 



r n (x;a,/3,S, N) 



( n) „=(*_*+„) vo U " W W (* + 1/ + * " J\0,-x+*-„ (2y + 5 - n + l) x - y * (gg) 
(-a + 5 - n) y {—P - n) y 



Moreover, 



(a + l)y((3 + S+l)y 

T^Nr n {x; a, (3, 5, N) = n(a + (3 + n + l)r n (x; a, (3, 5, N), 

{r n {-;a,P,S,N),r m {-;a,P,S,N)) VN = d„, m N — — . 

( n )N\{5 + /3 + l)„(a + 1)„ 

and r n (-; a, j3 1 5, N), n = 0, 1, 2, . . . , N, is a complete set of eigenf unctions for the operator T>n . 
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Proof. Just note that r n (x; a, ft, n, 8) is the unique function in V n satisfying the first order difference 
equation £ n tp = 0, with the initial condition ip(0) = 1; indeed, £ n tp — is equivalent to 

(x — a + 8 — n)(x — ft — n)ip(x) — (x + a + l)(x + (3 + 8 + l)<p(x + 1), x = 0, 1, . . . , n — 1. 

To get the norm, first note that 

a t mi2 W + l)n(« - S + 1)„ / *-n,($ -n)/2+l,-(0 + n), -(or- .5 + n),-n \ 
»i(-*)n \ (c5-n)/2,/3 + <5 + l,a + l,<5+l ' J 

+ l) n (5 -a- n) n (a + (3 + n+ 1)„ 



n!(<5 + /3 + l) n (a + l)„ 

where the last equality follows from ([55)1 (setting a = 5 — n, c = —(J) + n) and d = — (a — 6 + n)). Then 
using flSJ), (|30f and (f35|) , one can get easily the expression for the norm of r n (-;a, f3,8,N). 

□ 

In order to give the classical 4-F3 hypergeometric expression for r n , we will use the following identity 
of Whipple (see [TJ, Theorem 3.4.4): 



7^6 



a, a/2 + 1, b, c, <i, e, — n \ 
a/2, a — 6 + 1, a — c + 1, a — d + 1, a — e + 1, a + n + 1 ' / 

= ( a + i) n ( a -d-e+l) n ^ / o-6-c+l,d,e,-n \ 

(a-d+l)„(o-e + l) n 43 y a -6+l,a-c+l,d + e-a-?i ;1 / ' 
Theorem 7.3. The Racah polynomials have the following 4F3 hypergeometric expression: 

( a + P + n + l,-x,x + 6 - N,-n \ 
r n (x;a,/3,8,N) = 4F3 -i • 

V + 8 + l,a + l,-N ' / 

Proof. Consider the following elementary transformation formulas: 

(s + y + 5 - N) y - x+N - n (2y + 8 - n + l) x - v = ± — — ' — — -j- — (»" 

(x + o — N) y {6 - n){{o -n)/2 + lj y 

(n-y)(y) _ (N-X-n + l) n {-x) y {-n) y (8 + l) x 



(y + 8+l) a 



O " (N-n + l) n (N-x-n + l) y yV (8+l) y ' 

(x + 8-N) N+1 



(x + S- N) y - x+N - n = (5 - n) y 



(8 - n) n 



+X + 1 



(the last formula is valid, and must be proved, both for x > N — n and for x < N — n). Applying those 
formulas to (|36|) . we get easily the following >jFq expression for the Racah polynomials: 

, q r. {N -x-n + l) n (x + 8-n + l) n 

r n (x;a,f3,8,N) = ^ — ,',} x ~t\ -* 

(N — n+ l) n (8 -n + 1)„ 

(8 — n, (8 — n)/2 + 1, —f3 — n, —a + 8 — n, — x, x + 8 — N, —n \ 
1 • 
(8 -n)/2,S + + l,a + l,8-n + x + l,N -n-x + l,8 + l ' / 

Then Whipple's identity ([57]) gives immediately the 4F3 expression. 

□ 
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Remark 7.4. We have treated the case 7 = — N — 1 of the theory of Racah polynomials. The standard 
notation is r n {x; a, j3, 7, 6). 
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